Condensation in a zero range process on weighted scale-free networks 
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We study the condensation phenomenon in a zero range process on weighted scale-free networks in order 
to show how the weighted transport influences the particle condensation. Instead of the approach of grand 
canonical ensemble which is generally used in a zero range process, we introduce an alternate approach of the 
mean field equations to study the dynamics of particle transport. We find that the condensation on scale-free 
network is easier to occur in the case of weighted transport than in the case of weight-free. In the weighted 
transport, especially, a dynamical condensation is even possible for the case of no interaction among particles, 
which is impossible in the case of weight-free. 
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I. INTRODUCTION 

Condensation, a concept originally introduced by Bose and 
Einstein to explain the particle condensation in momentum 
space, is an intriguing phenomenon observed in real space, 
such as jamming in traffic communication , bunching of 
buses 031], and many mass transport models |4fl. In conden- 
sation, a finite fraction of particles may be condensed onto 
a single site. Nowadays, such phenomenon is found in the 
structure of complex networks, which characterize many nat- 
ural and manmade systems, such as the Internet, airline trans- 
port system, power grid infrastructures, and the world wide 
web ||5lla]. For example, Bianconi and Barabasi mapped the 
fitness model to a Bose gas by assigning an energy to each 
node and found that the fittest node can attract a finite frac- 
tion of all links [7, 8]. Besides the static condensation of the 
links of complex networks, the dynamical condensation on 
scale-free (SF) networks is found recently that the particles 
may completely condensed on the hub in the transport of zero 
range process (ZRP) where the interaction only occurs when 
the particles are stay at the same node 1191 1 1 Qfl . 

The previous studies of ZRP are mainly focused on regu- 
lar lattice 0, [U [H]. It is found that, in the steady state of 
one dimensional lattice, a finite fraction of the total mass will 
condense onto a single site when the global mass density is in- 
creased beyond a critical value. However, most of the realistic 
networks are not regular but SF Comparing to the reg- 

ular lattice, SF network is heterogeneous with power-law de- 
gree distribution P(k) ~ fc -7 and can be characterized by the 
Barabasi-Albert (BA) model with j = 3 lf]~3ll and its related 
modified models (H [H [M [HI [Hf] . It is pointed out that the 
structure inhomogeneity is an important factor in understand- 
ing many dynamical processes in SF networks lfl9i [20L Eltl . 
such as rumor propagation, virus spreading, and searching in- 
formation etc. For studying the influence of structure inho- 
mogeneity on ZRP, Noh et al. considered ZRP in BA model 
and found that the inhomogeneous structure makes particles 
to condensed on the hub when the jumping rate S is smaller 
than the critical point S c [HI [T^]. For a specific particle, its 
behavior can be considered as a random walk; and for a large 
number of particles, their behaviors are a kind of diffusion 
processes. As the diffusion makes the number of particles on 



a node fluctuate, the usually used approach to deal with the 
transport statistical properties in ZRP is the grand canonical 
ensemble, which can easily give the mean particles on each 
node. Recently, Ref. [4] uses the approach of canonical en- 
semble to explore the condensed phase and analyze the mass 
distribution. Here we introduce an alternate approach to deal 
with the ZRP, i.e., the mean field equations, and find that our 
approach will give the same results as that given by the ap- 
proach of grand canonical ensemble |9jll0[] when it is used to 
the same situations. Moreover, we use this approach to study 
the case of weighted transport and find that the added weight 
may push forward the critical value S c to a larger value 6' c , and 
S' c may be even over unity! 

Let's call the network with the same weight to transfer par- 
ticles on each link as weight-free network. Besides weight- 
free network, there is also weighted network where each 
link has different weight to transfer particles, which yields 
a more realistic description of communication in real net- 
works. Refs. [91 LL0|] investigated how the network structure 
influences particle dynamics. Our motivation here is to study 
how the weight of link influences the particle dynamics. In 
the weighted network, each link or node is associated with a 
weight Wi. The weight may represent the intimacy between 
individuals in social networks, or the bandwidths of routers 
and optical cables in the Internet. The strength of a node, Sj, 
is defined as the sum of all the weights of links of the node. 
It is found that the strength is strongly rely on its degree with 
s(k) ~ k a ll22l l23i 124(1 where a is different constants for 
different networks. For example, a is unity for the science 
cooperation network where the weight of a link between two 
scientists is given roughly by the frequency of their collabo- 
ration and 1.5 for the world-wide airport network where the 
weight is taken as the total number of passengers of the direct 
flights between two connected cities. The previous results on 
condensation focus on the case of a = In this paper, 

we consider the particle transport on the weighted network, 
i.e., the case of a > 0. Our results show that the condensation 
is easier to occur in the case of a > than that of a = 0. 
Especially, the condensation may occur even when the inter- 
action among particles does not exist, which is impossible in 
the case of a — 0. 

The paper is organized as follows. In Sec. II, we briefly 



2 



review the approach of grand canonical ensemble and its re- 
sults of condensation on SF networks with weight-free. Then 
in Sec. Ill, we give the dynamical mean field equations for 
the weighted transport on SF networks. Its limit will give the 
results of weight-free transport. In Sec. IV, we make numer- 
ical simulations to confirm the predictions given in Sec. III. 
Finally, the conclusions are given in Sec. V. 



II. THE APPROACH OF GRAND CANONICAL 
ENSEMBLE FOR THE CASE OF WEIGHT-FREE 
TRANSPORT 

The approach of grand canonical ensemble was first applied 
in the ZRP of one-dimensional syste m fl [Til [l2l l25ll and 
recently used in the SF networks | p, flQfl . Here we briefly 
review the results obtained in J^l Tioll . Suppose N parti- 
cles are randomly put on a network of L nodes and each 
node i can be occupied by any integral number of particles 
rii = 0, 1, 2, • • ■ ,N. Because of the interaction among the 
particles at the same node, some of the particles will jump out 
of the node and hops to other nodes, making the particle redis- 
tribution among n\, ria, • • ■ , n.£. Hence a microscopic config- 
uration is represented by n = m, n^, ■ ■ ■ ,til. The particles at 
node i with > may jump out with the jumping rates p(rii) 
and hop from node i to one of its neighbors j along the link 
with the hopping probability Tj<_i, In the SF networks with 
weight-free, Tj<_j is taken as 1/ki if i and j are linked and 
otherwise, i.e., a particle jumping out of the node i is allowed 
to hop to one of its neighboring nodes selected randomly, and 
the jumping rate is taken as p(n,) = nf. 8 = means only 
one of rii particles will jump out per each time, indicating the 
particles are attracting each other. In the case of 8 = 1, all the 
rii particles will jump out, implying they are moving indepen- 
dently and the system reduces to a noninteracting system of 
N random walkers. 5 > 1 means there is repulsive interaction 
among particles of one node, while 8 < 1 means the attractive 
interaction among particles. With enough time of evolution, 
the system will reach a stationary state. The mean occupation 
number in the stationary state is 



_ dlnFs(x) 

TfliyZ) — X ^ \x—zki 



(i) 



where z denotes the fugacity and 



71=0 v ; 



(2) 



It is found that the complete condensation occurs at 6 = 0, 
where the whole fraction of particles is concentrated at the 
hub |9, 10]. While there is no condensation at 5 = 1. For the 
case of 8 > 0, there is a critical 



Sc. = 1/(7 - 1) 



where the condensation will occur for 8 < 8 C and does not 
occur when 8 > 8 C . There is a crossover degree k c for 
< 8 < 8 C , which is defined as the degree with the aver- 
age occupation number mi = 1. rrii will be greater than one 



for the nodes with ki > k c and smaller than one for the nodes 
with ki < k c . There will be a condensation when there exists 
a finite k c > and no condensation otherwise. For i5 < 8 C , k c 
can be expressed as 



(Inkjyiax^ 



( k r , 



y-6/s c 



for 8 = 8 C 
for 8 < 8 C 



(4) 



where the nodes with k > k c will be the core of condensation, 
i.e., most of particles will hop to those nodes with k > k c . 
Furthermore, there is a scaling 



rrii = Gs(ki/k c 



(5) 



where the scaling functions behaves as Gs(y << 1) ~ y and 

G s ( y )(y > 1) ~ y 1/s - 



DYNAMICAL MEAN FIELD EQUATIONS FOR THE 
CASE OF WEIGHTED TRANSPORT 



III. 



According to the strength distribution s(k) ~ k a 11221 12311 . 
the strength is a nonlinear function of k and the nodes with 
the same k have the same strength. The strength will be dif- 
ferent for the ki neighbors of node i. When a jumping out 
particle hops to one of its neighbors, it will not choose all the 
links of the node on equal footing but with different rates to 
different neighbors. Therefore, we here assume the hopping 
rate Tj^i ~ fc". Making renormalization, we have 



TV 



k'i 



es 4 K r 



(6) 



where Bi denotes the set of neighbors of node i and Eq. (O 
will return to the case of weight-free when a = 0. 

We still suppose here that the jumping rates p(rii) satisfy 



and restrict 8 to < 8 < 1. With the evolution, 



rii will change with time and may be different for different 
nodes with the same k. For using the mean field approach, 
we transfer the description of n, for each node to the descrip- 
tion of the mean occupation number (t) for the nodes with 
the same k, i.e., rrik{t) is the average of all the rii of the nodes 
with degree k. Hence mu (t) is no longer necessary to be an in- 
teger. Correspondingly, we transfer the jumping rate p(rii) to 



After this transformation, let's see how many 



particles of a node can jump out at each time step. When 
mfc < 1, m s k > rrik and the jumping rate is in fact mk but not 
m s k as we don't have so many average particles mt at those 
nodes with degree k. When m& > 1, m s k < and hence the 
jumping rate is m 5 k . Except the aspect of jumping out, at the 
same time, a node accepts particles from its neighbors. The in- 
coming particles can be classified into two parts: One from the 
node with rrik < 1 and the other from the nodes with > 1. 
The incoming particles from one neighboring node with de- 
( 3 ) gree k' is P{k'\k)m k , {t)k a / £\ B fc" when m k ,(t) < 1 



and P{k'\k)m 5 k ,(t)k a / 'Y,j&B , k f when m k'{t) > 1, where 
P(k'\k) is the conditional probability for a link which rep- 
resents the possibility for a node with degree k to connect a 
node with degree k'. We use ko and k max to represent the 
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minimum and maximum degree of the network, respectively, 
and use k c (t) to denote the degree for m k (t) = 1. In the BA 
model, ko is a constant and k max ~ L 13 with /3 = 1/(7 — 1). 
By all these quantities, we have the mean field equations for 
the evolution of m k (t) 



dm k (t) 
dt 



fee 



-TOfe 



(t) + k(J2p(k'\k)m k ,(t) 



E 



J2p(k'\k)mi,(t) 



E 



h 



-),m k (t) < 1 



= -m{{t) + k{f^p{k'\k)m k ,{t) 



-),m k {t) >1(7) 



h 



where Ejes / ^j* denotes the sum to all the neighbors of the 
node with degree k'. Obviously, m k (t) depends on p(k'\k). 
For getting the detailed form of p(k'\k), here we consider 
the BA model as the underlying network. As the BA model 
is non-assortative mixing, its conditional probability satisfies 
p{k'\k) — k'p(k')/ < k > 11261 l27il . Borrowing the definition 
of the average degree of the nearest neighbors of a node with 
degree fc, i.e., k nn (k) = J2 k < k'p(k'\k), we have the average 



'Kn(k) = / P (k'\k)k' a dk' 
Jk 

= <k a+1 



> I < k > 



(8) 



It is easy to see that k" n does not depend on the degree k. 
Hence, we have Ejes fc , k j = ^'C an( l Eq. becomes 



dm k (t) 

dt 
dm k (t) 

dt 



= -m k (t) + k a+1 A{t), m k {t) < 1 

= -rn s k (t) + k a+1 A(t), m k (t) > 1 (9) 



where A(t) = Ejfc p(k')m k ,(t) + (*))/« 

Eq. © is the evolution equation of particle and will reach 
the stationary state when it evolves enough time. In the sta- 
tionary state, we may solve Eq. (O by dm k (t) /dt = to get 
the stabilized m k . By doing this we have 



m k 



k a+1 A, m k < 1 



m k = (k a+1 A)^ 5 , m k >l. 



(10) 



where A = )m k , + p{k') m {,)/{< k > 

k" n ) does not depend on the time t. We have A = k c ' a+1 ' 
from the condition m k = 1 for k = k c . Instituting A into Eq. 
( [Tol l we have 



m k = (k/k c ) a+1 , m k < 1 
m k = (k/k c Y a+1 ^ 5 , m k > 1. 



(11) 



Obviously, Eq. ( fTTT ) will go back to Eq. (O when a = 0. 

Recall that in the case of one dimensional lattice, one 
necessary condition for the condensation is that its density 
p = N/ L should be larger than a critical value p c > 0. While 
in the BA model of weight-free, there is no p c for p. The 
condensation occurs at any finite value of the particle density, 
i.e., p c — 0, and the particles can completely condensed at the 
hub or a few high-degree nodes ||9l [l0[]. Here the limit case 
of a — is the BA model of weight-free, hence we expect 
p c = for a > 0, which will be demonstrated later. For con- 
firming it, in the case of weighted transport we set the particle 
density p < 1. In the limit of L, N — > 00 we have m k < 1 
for all nodes if there is no condensation; otherwise, part of m k 
will become greater than unity. Therefore, the k c for m k = 1 
is nothing but the crossover degree. The number of condensed 
particles is 



Nc, 



N 



(k/k c ) 



(a+l)/S 



p(k)dk 



(12) 



In the condensed state, the nodes with k c < k < k max 
have the capacity to take an infinite number of particles. The 
non-divergence of the integration in Eq. (TTZb gives 



Si 



a + 1 



(13) 



The condensation will occur for the case of S < 5' c . From Eq. 
([T3l one can see that in the case of a = 0, 8' c = 1/(7 — 1) 
returns to the case of weight-free of Eq. ©. By N con /N ~ 
O(l), we have 

r k max 

k (a+l)/S _ / k^ a+1 ^ 6 k-^dk, (14) 

Jk c 



and further we have 



(lnk m ax) Sc , for 6 = 5'c 
(hna,) 1 -^^, for 5<5' c ( ' 

Obviously, it has the same form with Eq. (0]i except S c is 
replaced by 8' . From Eq. (fT~5T > one can see that k c does not 
depend on the particle density p, indicating the condensation 
can occur at any p. Therefore, we also have p c = for the 
weighted transport. 

For the situation of 5 < 8' c , from Eq. ( fTTT i the particles 
accumulated at the hub is 



m h ub ~ (k ill (I X 

- (k S/S '<')( a + 1 V S - T 



(16) 



indicating rrihub scales linearly with L and hence there is con- 
densation at the hub. And for the situation of 8 > 8' c , the 
particles at the hub is 



mhub 



k (a+l)/S _ L S' c /5 



(17) 



indicating rrihub scales sublinearly with L and hence there is 
no condensation. When there is condensation, from Eq. ( ITTb 
we have 



m k r~j k v 



(18) 
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with i] = (a + l)/6 for 8 < S' c and mu > 1. The scaling T) 
will decrease to (a + l)/5' c = 7 — 1 when 8 increase to 8' c , 
which gives a critical r\ c = 2 for 7 = 3. Therefore, we can 
judge the occurrence of condensation by checking the scaling 
?y for mk > 1, i.e., there is condensation if r\ > 7 — 1 and no 
condensation otherwise. 

Moreover, from Eq. ([T3T > we can see that, for BA model 
with 7 = 3, it is possible to have S' c > 1 for a > 1. In 
this case the condensation will occur for the whole range of 
< 8 < 1, which is impossible in the case of weight-free. 
For example, taking 8=1 for the case of weight-free, all the 
particles at node i will jump out of the node per each time, 
resulting a uniform distribution of particles and hence no con- 
densation. However, in the case of weighted transport, the 
situation is totally changed. During the evolution, all the par- 
ticles will jump out at each time and hop with the favorite to 
one of the neighbors with the largest degree. Gradually, most 
of the particles will move to the larger and larger nodes, and 
finally to the hub around. After that, all the particles at the 
hub will jump out per each time when 5 = 1, but at the same 
time, most of the particles of the hub's neighboring nodes will 
hop to the hub, resulting in a dynamical equilibrium. Namely, 
6 = 1 for 6' c > 1 is a kind of completely exchange conden- 
sation, which is different from the case of 6 < 1 where most 
of the previous particles at the hub will still stay at the hub. 
This kind of complete exchange condensation also exists for 
the case of 1 < 6 < 8' c where the jumping rate p(mk) = mj. 
is the same with that of 5 = 1. 



IV. NUMERICAL SIMULATIONS 

In numerical simulations, we first construct a BA grow- 
ing network with the size L = 2000, the average degree 
< k >= 6, and the degree distribution P(k) ~ fc~ 3 accord- 
ing to the algorithm given in Ref. d 1 1 811 - Then we put 
N = 1000 particles at the L nodes randomly, i.e., the particle 
density p = 0.5. Hence m^it = 0) is uniform for different k, 
see the horizontal line for the case of t = in Fig. [TJ At each 
time step, we let the N particles evolve according to the jump- 
ing rate p(rii) and the hopping rate 1^<_$ in Eq. ([6}. Namely, 
we choose min(n,i, nf) particles from the n, particles of node 
i and let them hop to the neighbors of node i according to Eq. 

where min(m, nf) means taking the smaller one from m 
and nf. By this way, we find that the numerical simulations 
completely confirm the evolution Eq. (0 and the stationary 
distribution Eq. (fTTT l. For example, from Eq. (0 we have 
n%k(t = 1) = A(0)k a+1 ~ k a+1 , indicating a power-law re- 
lation between and k for all the k at t = 1, Our numerical 
simulation has confirmed it, see the straight line for the case 
of t = 1 in Fig. [TJwhere the parameters are taken as a = 0.4 
and 8 = 0.2 < 5' c = 0.7. This algebraic relation will be 
kept until the emergence of some nik(t) > 1. After that, we 
know from Eq. (0 that the nodes with m^it) < 1 for small 
k will keep the algebraic relation mfc(i) ~ A(t — l)k a+1 but 
the nodes with m,k(t) > 1 for large k will depends on both 
k and mu(t — 1) because of the nonlinearity of the first term 
m s k (t) in Eq. (0. The numerical confirmation is given by the 




FIG. 1 : Evolution of particle distribution in weighted transport for an 
arbitrary a = 0.4 and S = 0.2 < S c = 0.7. 



curves for the cases of t = 8 and 64 in Fig. Q] It is easy to 
see from these two curves that the part with smaller k has the 
same slope with that of the curve of t = 1 and the part with 
larger k is complicated. After further evolution, the system 
will reach its stationary state. The lowest lines in Fig. Q]show 
the numerical results for t = 1024 and 40000 respectively. 
Obviously, the case of t = 1024 is overlapped with the case 
of t = 40000, indicating the stationary state is reached. 

In the stationary state, from Eq. (fTTl i we know that the 
stabilized mjt will have two slopes and both depends on the 
weight scaling a, i.e., slope a + 1 for < 1 and slope 
(a + l)/5 for mfc > 1. The numerical simulations have com- 
pletely confirmed this. Fig. |2ja) shows the results for a = 0.4 
and S = 0.2, 0.5, and 0.8 respectively. By Eq. < TT~3T > we get 
6' c = 0.7, so the case of 8 = 0.2 and 0.5 should have conden- 
sation and the case of 8 = 0.8 > 8' c no condensation. This can 
be confirmed by checking the scaling r\ in Eq. (TT~8T > . It is easy 
to see from Fig. |2J a ) mat me slopes are 7.0 and 2.8 for the 
cases of 8 = 0.2 and 0.5, which are greater than t] c = 2, and 
1.75 for the case of 8 = 0.8, which is less than 2. These val- 
ues also confirm the prediction of Eq. ([Til in which we have 
the slopes a + 1 = 1.4, and (a + l)/8 = 7, 2.8, and 1.75 for 
8 = 0.2, 0.5, and 8 = 0.8, respectively, see the solid reference 
lines with different slopes s in Fig. |2fa). On the other hand, 
from Eq. (fT~5b we have lnk c ~ C8. The crossover degree k c 
can be measured from the intersection between the dotted line 
TOfc = 1 and the measured curves in Fig. |3Ja). By this way 
we can get k c for different 8. Fig. |2jb) shows the result. It is 
easy to see that Fig. |2jb) is a straight line with negative slope, 
which qualitatively explain Eq. ( fT3T >. 

By the fact that 77 > r] c means condensation, we may deter- 
mine S' c numerically. For a given a, we increase 8 gradually 
and check 77 for each 8. The critical 8' c can be identified as 
the value that separates the 77 > 2 and 77 < 2. By this way, 
we obtain 8' c for different a. Fig. [3] shows the result. Ob- 
viously, it is a straight line with slope 0.5 and thus confirms 
Eq. (IT3T >. When a is over unity, however, we cannot obtain 8' c 
numerically as 8 is limited to unity for the jumping rate p(n). 

Furthermore, from Eq. (fTTT l we can see that the two slopes 
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FIG. 2: (a) Stabilized rrik versus k for a — 0.4 and S' c = 0.7 where 
the lines with "squares", "stars", and "circles" denote 5 — 0.2, 0.5, 
and 0.8, respectively. The drawn lines for slope s = 1.4, 7.0, 2.8, 
and 1.75 are for reference, (b) Crossover degree k c versus a. 




FIG. 3: The critical parameter S' c of condensation versus the weight 
parameter a. 

for nik < 1 and m/- > 1 will become the same when 5=1. 
This phenomenon is also confirmed by our numerical sim- 
ulation. Fig. [4] shows the result where the measured slope 
2.37 (the solid line) is consistent with the theoretical predic- 
tion a + 1 = (a + l)/6 = 2.4 and the fact that 2.37 > 2 
confirms the condensation for 6 = 1 < S' c . 



V. DISCUSSIONS AND CONCLUSIONS 



The condensation in ZRP was usually discussed in mass 
transport and only a little attention had been paid to the inter- 




FIG. 4: The stabilized vtik versus k for ot = 1.4 and 8 = 1 < 8' c = 
1.2. 



acting dynamical systems on SF networks. The work of Noh 
et al. makes ZRP on SF networks become reality. As com- 
plex network has very rich features, Noh's work is very im- 
portant but only the beginning of ZRP on complex network. 
There are a plenty of factors in networks which may affect the 
condensation, such as the distribution, clustering coefficient, 
assortativity, and weight etc. This paper investigates one of 
the factors. Our result show that the weight may significantly 
influence the condensation and make it possible for the con- 
densation to occur even for S — 1 when 5' c > 1, which is 
impossible in the weight-free network. We are continue along 
this way. 

In conclusions, we have discussed the dynamical conden- 
sation on SF network with weighted transport. We have in- 
troduced a set of dynamical mean-field equations in ZRP to 
describe the evolution of particle and find that the particle con- 
densation depends on not only the jumping rate, but also the 
weight of transport. Our results show that the mean field ap- 
proach is completely equivalent to the grand canonical ensem- 
ble in the situation of a = and the condensation is easier to 
occur in the weighted network than that in the network with 
weight-free. The critical value of condensation in weighted 
network is S' c = (a + l)/(7 — 1), which is larger than the 
value 5 G = 1/(7 — 1) of weight-free. 

This work was supported by the NNSF of China under 
GrantNo. 10475027, by the PPS under Grant No. 05PJ14036, 
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